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Abstract. Large scale numerical simulations are used to study the elastic dynamics 
of two-dimensional vortex lattices driven on a disordered medium in the case of weak 
disorder. We investigate the so-called elastic depinning transition by decreasing the 
driving force from the elastic dynamical regime to the state pinned by the quenched 
disorder. Similarly to the plastic depinning transition, we find results compatible with 
a second order phase transition, although both depinning transitions are very different 
from many viewpoints. We evaluate three critical exponents of the elastic depinning 
transition. (3 — 0.29 ±0.03 is found for the velocity exponent at zero temperature, and 
from the velocity-temperature curves we extract the critical exponent 6^^ = 0.28±0.05. 
Furthermore, in contrast with charge-density waves, a finite-size scaling analysis 
suggests the existence of a unique diverging length at the depinning threshold with 
an exponent v = 1.04 ± 0.04, which controls the critical force distribution, the finite- 
size crossover force distribution and the intrinsic correlation length. Finally, a scaling 
relation is found between velocity and temperature with the /3 and 5 critical exponents 
both independent with regard to pinning strength and disorder realizations. 
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1. Introduction 

The physics of disordered elastic systems is relevant for numerous systems, e.g. interfaces 
hke magnetic [1] or ferroelectric [2] domain walls, contact lines [3], crack propagation 
[4], or periodic structures like charge- density waves [5] (CDWs), vortex lattices in type 
II superconductors [7j, colloids [lO] or Wigner crystals [H]. It is essential to understand 
the response of these systems to an external driving force, such as current-induced 
Lorentz force for vortices, electric field for charge- density waves, or magnetic field for 
domain walls. The competition between elasticity and disorder in all these systems 
leads to a great variety of phases. In particular the system remains pinned up to a 
critical driving force and starts sliding above F^. The transition from a pinned state 
to a sliding one is known as the depinning transition. It has been suggested by Fisher 
|12] . in the context of CDWs, that the depinning transition could be regarded as a 
standard equilibrium critical phenomenon where the reduced force f = {F — F^) / F^. and 
the velocity v would respectively act as the control parameter and the order parameter. 
This idea has proven useful in many other systems. This analogy suggests the existence 
of critical exponents, in particular /3 which describes the power-law response of the 
velocity at zero temperature v ~ . However the analogy with standard critical 
phenomena is limited since the depinning transition is by nature a non-equilibrium 
transition. For example, unlike equilibrium phase transitions no divergent steady-state 
correlation length scale exists for an elastic line when approaching the depinning from 
below [13] . Connections between the depinning transition of elastic interfaces in random 
media and non-equilibrium phase transitions into absorbing states are reported in the 
literature [15]. The pinned phase where the dynamics is frozen may be considered as 
an absorbing state in which the system may enters but cannot leave. Two universality 
classes have been proposed for the interface growth. The linear growth of the interface 
is described in the quenched Edwards- Wilkinson equation within the random field 
Ising class [191 1201 [21], while Kardar-Parisi- Zhang nonlinearities in this equation make 
the depinning transition closely related to directed percolation [22]. Non-equilibrium 
phase transition into absorbing states has also been shown in periodic systems with 
quenched disorder like colloids: when periodically sheared a transition was observed 
from a diffusive liquid to a pinned state [21], while in the case of driven colloids the 
transition separates plastic flow from a pinned state [25]. The critical exponents found 
in both studies are consistent with aborbing state transitions in the universality classes 
of directed percolation or conserved directed percolation. 

When the strength of the underlying disorder is strong, the theoretical description 
of the phenomenon is difficult and the question about the nature of the depinning tran- 
sition remains an open problem. In periodic systems with a 2D displacement field of 
dimension (A^ = 2) (e.g., superconductor vortices, colloids and Wigner crystals), simu- 
lations [26| [M| [37] have shown that this situation leads to dislocations and plasticity: at 
the depinning threshold tearing of the lattice appears where regions of pinned particles 
coexist with particles flowing around them. When the pinning of the underlying sub- 
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strate is weak the statics and dynamics are dominated by the system elasticity which 
favors order. After the work of Fisher, further analytical and numerical works have 
refined the analogy with equilibrium critical phenomena to define universality classes 
and critical exponents (see e.g. [38] and references therein). However in numerical sim- 
ulations of periodic systems in ti = 2 various exponents have been found. /3 ~ 2/3 is 
found for CDWs [3^, Wigner crystals [IQ], and colloids [37], while other works found 
/3 ^ 0.5 and (3 = 0.92 ± 0.01 jl2] in colloids, and (3 = 0.35 [13] in stripe systems. 
For superconductor vortices in d = 3 the value /3 = 0.65 is found [B], while in d = 2 a 
surprising exponent (3 > 1 has been reported [15]. 

In this paper, we show a detailed study of the critical behavior of superconductor 
vortices above the elastic depinning transition, including a finite-size scaling analysis 
and the determination of three critical exponents. The elastic depinning transition 
investigated in this paper is very different from the plastic depinning transition studied 
in a previous paper [16]. In the plastic case tearing of the lattice at the depinning 
threshold is observed with moving vortices flowing along changing interconnected 
channels around pinned regions of vortices. On the contrary, at the elastic depinning 
threshold, all vortices depin together and flow along well-defined rough static channels 
in a topologically ordered structure. Furthermore, the velocity is periodic in time at the 
elastic threshold, while it is shown to be chaotic at the plastic threshold [3l] with large 
broad-band noise at low frequency. The curvature of the velocity-force curve is also 
very different since the elastic depinning is caracterized by a velocity exponent /3 < 1 
whereas /3 > 1 for the plastic depinning. 

We perform large scale molecular dynamics simulations of 2D superconductor vortex 
lattices driven over weak random disorder. Our model belongs to the category of 
2D periodic systems with N = 2 and short-range interactions. It could model 3D 
superconductors (either conventional or layered) in an effective 2D regime, i.e., when 
the vortex line tension is high enough for the lines to remain straight. The numerical 
model is developed in section |2l The depinning transition at zero temperature is detailed 
in section |3] for different lattice sizes and disorder realizations. Above the depinning 
threshold a power law response vf>Fc ~ is found with (3 = 0.29 ± 0.03, and the 
transition seems to be governed by the correlation length divergence with an exponent 
u = 1.04 ±0.04. The same exponent governs the finite-size effects. In section H] we study 
the depinning transition at finite temperature. A power law response for the velocity 
vf=Fc ~ T^^^ is found and the critical exponent = 0.28 ± 0.05 is determined. A 
scaling relation for the driving force and the temperature is established. Varying system 
size, disorder realization and pinning strength, the scaling function and both /3 and S 
critical exponents show some degree of universality. The critical exponent values of our 
study are discussed in section 5. 
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2. Numerical model 

We study A'^, Abrikosov vortices interacting with Np random pins in the (x, y) plane. 
Periodic boundary conditions are used in both x and y directions. We consider the 
London limit \l ^s, where Al is the penetration length and is the superconducting 
coherence length, i.e. we treat vortices as point particles. As in |17], the LAMMPS 
classical molecular dynamics code [H] is used to integrate the newtonian equations of 
motion with a velocity Verlet algorithm |19]. We use this parallel code in order to 
simulate large system sizes, which is crucial to study the elastic depinning since the 
Larkin domains are larger than in the case of plastic dynamics. However the LAMMPS 
code includes an inertial term proportionnal to the particle mass. The existence of 
a vortex mass has been discussed in the literature [50| 155] . In general vortices are 
considered as massless objects although the vortex mass may be significant in the clean 
limit. In recent experiments |55], it is shown that Nb vortices close to Tc ~ 9K have 
an inertial mass which is in the intermediate range between the dirty and clean limit, 
showing that the inertial vortex mass should be a meaningful concept. 

We describe the time evolution of a vortex z at a position rj with the equation 

+ + Ff (t) 

where Vij is the distance between the vortices i and j located at and r^, rip is the 
distance between the vortex i and the pinning site located at r^, and Vj is the 2D 
gradient operator acting on r^. F^ = Fx is the Lorentz driving force in the x direction 
induced by an applied current, and rj is the viscosity coefficient. The Langevin force 
F**^ describes the coupling with a heat bath and can be expressed by a thermal gaussian 
white noise with zero mean and variance 

= 2vkBT6^,,6,,6it - t') 

where X, ^ = x,y are the cartesian components of F*^ and ks is the Boltzmann constant. 
The vortex-vortex pairwise repulsive interaction is given by a modified Bessel function 

U^^{r,j) = a,Ko{r,,/XL) 

and the attractive pinning potential is given by 

U^p(^rip) = -a,e-('--/^-)' 

where Rp the radius of the pins, and a„ and ap are tunable parameters. We fix the 
strength of the vortex-vortex interaction by setting a^, = 2.83 10~'^ALeo where eg is a 
characteristic energy per unit length. The number of pins is set to A'^ = A'^;; and their 
radius is Rp = 0.22Xl. The positions of the pins are chosen from a uniform random 
distribution. Each random sampling will be referred to as one "sample" or one "disorder 
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realization". The average vortex distance is Qq = -^l? and the vortex- vortex interaction 
is dealt with using a neighbor list method with a cutoff radius = 6.5 Xl. Since we 
are mainly interested in the permanent regimes we choose rj/m = 0.1 (where m is the 
vortex mass) in the LAMMPS code such that the inertial term is small compared to the 
viscous term. Care has been taken to check that such ratio ensures that the long time 
behavior of the second order Newton's vortex dynamics is identical to the overdamped 
dynamics limit usually computed for the superconductor vortices. In particular we 
find identical vortex trajectories and identical values of the critical exponents using 
overdamped dynamics on smaller systems. 

We use a unit system in which ?7 = l,AL = l,eo = l and = 1- Depending of the 
relative strengths of the vortex-pin and vortex-vortex interaction, the dynamics can be 
either dominated by elasticity or disorder. Several values of the relative disorder strength 
ap/«„ have been investigated. In sect ions ISTTl 13. 2 1 and 14. II we choose the relative disorder 
strength ap/a^ ~ 5.10"^, weak enough to enforce elastic behavior (plasticity is found 
above a^ja^ ~ 0.01). In section I^y2] we investigate several pinning strengths in the elas- 
tic regime, from a^ja^ ~ 2.10^^ to a^ja^ ~ 7.10"^. Molecular dynamics simulations are 
performed for several system sizes, from A^^ = 270 up to = 20000 vortices, but keep- 
ing fixed vortex and pin densities. Moreover, various rectangular shaped basic cells of 
size {Lx,Ly) have been investigated, from the almost square (L^^Ly) = {5,6\/3/2)n\L 
with n = 3,8,20 up to very elongated strips {Lx,Ly) = (400, 20\/3/2)Al. The strip 
geometry elongated in the driving force direction allowed the study of the critical de- 
pinning properties for large system sizes. Care has been taken to check that the basic 
cell anisotropy induced by the strip geometry does not alter the critical properties, and 
in particular the value of the velocity critical exponent (3. 

3. Zero temperature 

For several disorder realizations we start from a perfect triangular lattice at high ve- 
locity and the driving force is slowly decreased down to the sample-dependent critical 
depinning force i^'^ampie j^gjQ^ which the system is permanently pinned. We compute 
the mean critical depinning force Fc where the overline indicates average over disorder 
realizations. We show in FigJT^ the evolution of Fc with respect to the relative disorder 
strength Op/a^ for = 8000 vortices in a basic cell of size (L^., Ly) = (400, 20-\/3/2)Ai. 
The rapid increase in the mean critical depinning force indicates a crossover from elastic 
dynamics dominated by elasticity to plastic dynamics dominated by disorder: such a 
sharp increase was for example found in a previous work on a similar system [37] . In Fig. 
[lb the typical trajectories of the vortices at the elastic depinning threshold are displayed 
for otpjav ~ 5.10"^. The vortices flow in elastically coupled rough static channels and 
the structure is topologically ordered: all vortices depin together with the same mean 
velocity, which means that all vortices keep the same neighbors as they move. The 
dynamics is jerky and the velocity of the center of mass is periodic in time where the 
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Figure 1. a) Transition from elastic to plastic dynamical regimes shown by 
the evolution of the mean critical depinning force -Fc with respect to the relative 
disorder strength ctpja^, for = 8000 vortices in a basic cell of size {Lx,Ly) = 
(400, 2G\/3/2)Al. b) Typical trajectories of vortices at the elastic depinning threshold 
obtained for — 12000 vortices and ap/a^ w 5.10"'^. A snapshot of the vortex 
positions (filled circles) at a given time is superimposed to their trajectories. For 
clarity only a small part of the basic cell of size iL^,Ly) = (100, 120y3/2)AL is shown. 



period corresponds to the time for each vortex to replace its preceding neighbor in the 
same channel 



3.1. Velocity-force response 

For each system size {Lx, Ly) we define the reduced velocity, 

v = < v-"^{t) > /Fr^^p'' 

and the reduced force f = {F — iTfampic^^^sampic -y^^j-^g^g < Vx^it) > is the time averaged 
longitudinal velocity of the center of mass of the vortices and i^^^ampie jg ^^iq critical force, 
both measured for a given disorder realization, and the overline is an average over all 
disorder realizations for a fixed value of /. This approach has been used for example in 
numerical simulations of elastic interfaces in a disordered medium [56] . 

We set the relative disorder strength ap/a^ ~ 5.10"^ and we plot in Fig. [2] 
the reduced velocity v with respect to the reduced force / for two different system 
sizes: {Lx, Ly) = (100, 50-\/3/2)Ai averaged over = 21 samples in Fig. [2^, and 
{Lx,Ly) = (400,20V3/2)Al avera ged over = 14 samples in Fig. [2b. 

Three regions appear in Fig. [2l Region I is the manifestation of the finite size 
effects in the system whose signature is the single-particle regime [16], where v ~ /^/^ as 
shown by the lines of slope 1/2. Possible hysteretic depinning may be measured in this 
region for few disorder realizations, with different values F"^ and F^°^^ of the threshold 
force when increasing or decreasing the force. However, the width F^^ — of the 
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Figure 2. Plot of the vortex reduced velocity v versus the reduced force / for 
different system sizes: a) {Lx,Ly) = (100, 50\/3/2)Al for N = 21 samples b) 
(L^jLy) = (400, 20-\/3/2)Ai for iV = 14 samples. The solid lines indicate a slope 
of 1/2 corresponding to the single-particle regime. 



hysteresis decreases when the system size increases, which therefore confirms that such 
hysteresis phenomenon is just a finite size effect. In region II a power law regime f ~ 
with /3 < 1 is measured, which we identify with the critical regime of the continuous 
depinning transition. Finally, in region III the system is far above the critical depinning 
threshold and approaches the asymptotic linear behavior f ~ / obtained in the absence 
of disorder. 



/5 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



20 40 60 80 100 120 140 

Ly 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 







.o 

















100 



200 300 



400 500 



Figure 3. Depinning exponent /3 extracted from a power law fit in region II 
of the disorder averaged velocity- force response a) versus Ly for = lOOA^. 
Averages are computed over TV = (9,7,4,21,4) samples for respectively Ly = 
(2,7,18,50,120)v^/2Ai, b) versus L^ for various transverse sizes from Ly = 18V3/2Al 
to Ly = 120V3/2Al. Averages are computed over TV = (44,47,39,12,27) samples for 
respectively L^ ~ (15,40, 100, 200, 400)Ai. The error bars correspond to the standard 
error associated to the value of /3 extracted from the fit. 



For each system size we compute the depinning exponent /3 from a power law fit in 
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the region II of the disorder averaged v—f curves (see Fig.H]). Fig. shows the evolution 
of P with respect to the transverse size Ly for a fixed longitudinal size = lOOA^. We 
see that the value of (3 for Ly > 18-\/3/2Al becomes independant of the transverse size 
Ly. In particular /3 does not depend on the basic cell anisotropy since we measure 
identical values in square basic cells (L^^Ly) = (100, 120y3/2)AL. In Fig. Eb we show 
the evolution of /3 with respect to the longitudinal size L^ for various transverse sizes 
Ly > 18V3/2Al. It can be seen that /3 has reached a constant value for L^ > lOOA^. 
Taking the mean value of (3 computed for both L^ > lOOA^ and Ly > 18-\/3/2Al, we 
obtain the result (3 = 0.29 ± 0.03. 

Note that taking the mean of individual values Pi for each disorder realization, i.e. 
without averaging the velocity over disorder, gives a similar value /3 = 0.27 ± 0.04. 

3.2. Finite-size scaling 

Finite-size systems used in numerical simulations lead to rounding and shifting effects 
in second-order phase transitions that can be analysed within the finite-size scaling 
theory [57] . In addition to providing a clearer picture of the thermodynamic limit 
of the transition from finite size observations, this approach allows to extract further 
information about diverging lengths in the system. From a general theorem for 
disordered systems [58] a finite-size scaling length can be defined from the statistical 
properties of a large number of finite-size samples. The finite-size scaling exponent 
Uys characterizing the divergence of such a length at the threshold should satisfy the 
inequality >2/d. 

We first examine the critical force distribution for each system size, and determine 
their width AFc(L^) for the relative disorder strength a^/a^ ~ 5.10"^. In Fig. we 
display AF^lL^) versus L^ which scales as [59] , 

with z/pg = 1.09 ± 0.07, in agreement with the inequality Ups > 2/(i. 

Similarly, for each system size we identify the appearence of finite size effects with 
the crossover force -Fo-oss between the critical regime and the single-particle regime (see 
regions I and II in Fig. [2]). We plot in Fig. |3)d the width AF^rossiLx) of the crossover 
force distribution, which behaves like 

with the same exponent z/pg = 1.09 ± 0.02. 

Furthermore, we expect the correlation length ^ to behave in the critical regime as 
^ ~ (F — F^)~'^ where v is the intrinsic correlation length exponent, and F^ is the 
critical force for the infinite system. At the crossover force -Fcross the correlation length 
becomes of the order of the system size L^, so that the crossover force -Fcross varies with 
the longitudinal size L^ as 

-Across c 
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Figure 4. a) Width of the critical force distribution AFc{Lx) versus L^.. A power- 
law fit to the data gives the solid line of slope — l/z^ps from which we extract 
I'ps = 1.09 ± 0.07; b) Width of the crossover force distribution AFcross(-^x) versus 
Lx. A power-law fit to the data gives the solid line of slope —1/v-ps from which we 
extract z^ps = 1.09 ± 0.02. The data are the same as those shown in Fig. [3J). The error 
bars indicate the standard error associated to the value of the width extracted from 
the distribution. 



Fig. [5] shows Fcross — versus where Fcross is the disorder averaged crossover force. 
We determine from the asymptotic behavior of the disorder averaged sample critical 
force Fc for large Lx (see the inset of Fig. [5]). 




Figure 5. Plot of Fcross-- versus for the same data as in Fig. ID A power-law fit 
gives ly = 1.04±0.04. Inset: Averaged sample critical force Fc versus Lx from which we 
extract the asymptotic value F^ for large L^- The error bars represent the standard 
error of our samples. The large error bars for the two largest system sizes in the main 
figure come from the statistical errors which become of the order of i^cross — F^ . 

From a power law fit shown in Fig. El we extract the value z/ = 1.04 ± 0.04 which 
is consistent with z/pg computed previously. Therefore, a diverging length above the 
depinning threshold with an exponent u = z/pg seems to appear. Our results suggest 
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that this length coincides with the intrinsic correlation length ^. We show in Fig. E]the 
direct computation of the correlation length ^ extracted from the velocity correlation 
function fit with the functional form < Vx{0)vx{x) >~ Cox~'^e~^^^ (see [60]). We find 
that both fit parameters Cq and k are almost constant when F — varies. For a given 
disorder realization and {L^^Ly) = (100, 18-\/3/2)Al, Fig. |6] shows the evolution of ^ as 
a function oi F — F^. 
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Figure 6. Plot of ^ versus F — F^ for a given disorder realization and for a system size 
{L^,Ly) = (100, 18V3/2)Al at ap/a„ « S.IO-^. The error bars indicate the standard 
error associated to the value of ^ extracted from the fit < Vx{0)v^{x) >^ Cox~'^e~^^^ . 
The black arrow indicates the value of Fcross — F^. 

As expected, when F — )■ -F^oss finite size effects actually appear when ^ ~ L^, with 
Lx = IOOAl. Furthermore the correlation length ^ diverges as ~ ~ F^)~'^ for 
F > -Fcross in the critical regime, and the power law fit exponent gives z/ = 1.2 ± 0.2 
which nicely agrees with Fig. El Therefore, our results strongly suggest that, unlike 
CDWs [591 EH |63] , the intrinsic correlation length controls both the critical force dis- 
tribution and the finite-size crossover force. 



We can now write a scaling relation between the velocity, the driving force and the 
system size L^, at T = as. 

This scaling relation is relevant only close enough to the depinning transition and at 
large enough L^- 

We take the disorder averaged response v{f) for different system sizes = 
(40, 100, 200, 400)Al and various transverse sizes from Ly = 18a/3/2Al to Ly = 
120\/3/2A2,, and for the relative disorder strength ap/a^ ~ 5.10"'^. In Fig. [7] we plot 
vLx^'^ versus fL]J'^ using the value of /3 = 0.29 ± 0.03 and u = 1.04 ± 0.04, retaining 
only the points of the region I and II. A good collapse of the data onto a single curve 
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Figure 7. Plot of vLx^^ versus JL^J'^ for different system sizes Lx = 
(40, 100,200,400)Al at ap/a^ « 5.10^^. The value of the /3 and v exponent is 
/3 = 0.29 ± 0.03 and v = 1.04 ± 0.04. 

is obtained. The behavior of the scahng function g{x) is g{x) ~ for a; ^ 1, and in 
contrast to ordinary phase transitions g{x) — )■ for x — )■ 0. This unusual property has 
aheady been observed in numerical simulations of the depinning transition of driven 
interfaces [64]. It is due to the single-particle regime which can partially hides usual 
finite-size effects. 




4. Finite temperatures 

4^.1. Critical exponent 5 

In this section we study the finite temperature effects on the depinning transition. For 
-F < Fc at T > the vortex mean velocity does not vanish since thermal fluctuations 
provide sufficient energy to overcome local energy barriers: this results in a rounded 
depinning transition. Using the analogy with standard critical phenomena, a power law 
response vp=p^ ~ T^/^ can be defined at F = F^. We show in Fig. [S] the typical reduced 
velocity-temperature response that we obtain. Two different behaviors can be seen in 
the f (T) response. Below the critical force, v goes to zero faster than a power law 
resulting in concave v{T) curves, while above the critical force v approches a nonzero 
limit as T goes to zero resulting in convex v{T) curves with an horizontal asymptote 
on the left. The effective critical force F* is defined as the force at which the convexity 
changes. In agreement with a second-order phase transition, we can extrapolate at F* 
a power law response Vp^j^ ~ T^^^ (dashed line in Fig. [S]) from which we measure the 
critical exponent = 0.28 ± 0.05, where the error bar is deduced from the different 
lines one can draw to extrapolate the power-law behavior. 
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We would like to emphasize that this new determination of the disorder averaged 
critical force is consistent with the value of Fc analysed at T = for the same samples 
in section [3l Going back to the v{f) curves at T = and using the effective critical 
force F*, we find (3 = 0.28 ± 0.05 which nicely agrees with the value /3 = 0.29 ± 0.03 
extracted from the zero temperature simulations of section [3l This larger uncertainty 
comes from the smaller number of samples simulated at finite temperatures. 

Finally, note that we observe an identical mean value = 0.28 when taking the 
mean of individual values obtained from the velocity-temperature curves for each 
disorder realization. 




Figure 8. Reduced velocity v versus temperature T, averaged over N = 6 samples 
for fixed values of / = (10,4,3,0,-2,-7) x 10~^. The system size is {Lx,Ly) = 
(400, 20-\/3/2)Al and the relative disorder strength is ap/av ~ 5.10"^. The dashed 
line is the extrapolation of v(T) at the effective critical force F* from which we extract 
the value of 6~^. The error bars correspond to the standard error. 



4-2. Different pinning strengths - Scaling law 

For different system sizes with > IOOAl and Ly > 18V3/2Al, we investigate sev- 
eral values of the relative pinning strength ap/a^ ~ (2,3,5,7) x 10~^. The number of 
samples is = 1 for ap/a^ ^ 2.10"^ and = lOOOA^, = 2 for ap/a^ ^ 3.10"^ 
and = 400Xl, iV = 8 for Op/a^ ^ 5.10-^ and = (100,400)Al, iV = 4 for 
oipjoiv ~ 7.10"^ and = 400Al. For each pinning strength, we compute both the crit- 
ical exponents (3 and and we display their evolution respectively in Fig. ^ and Fig. 
|9]3. The values of (3 and are found to be independent of the relative pinning strength. 

Because of the power law responses of the velocity with T and /, we assume that 
the velocity w is a generalized homogenous function of driving force / and temperature 
T. We want to express the relation between v, /, and T in terms of dimensionless 
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quantities and to make this relation independent of the pref actors in the two power law 
responses. We define Vq and Tq as 

Vf>0,T=0 = Vof'^,Vf=o = Vq f — j 

and we introduce the dimensionless velocity v = v/vq and the dimensionless temperature 
T = T /Tq. We can then define the following scaling ansatz, 

v\f\-^ = S^{f\f\-^') 

where S{x) is the scaling function with two branches 5*+ and corresponding to / > 
and / < 0, respectively. The power law dependencies 'y/>o,T=o = and Vf=o = T^/^ 
imply that S{x) has the following asymptotical behavior: 

lim5'+(x) = 1 and lim x~^^^ S±{x) = 1 

The S+ branch goes asymptotically to the horizontal axis for T — )■ defining a driv- 
ing dominated regime, while the S± branches have both an oblique asymptote with 
slope for / — ?■ defining a temperature dominated regime. The change of variables 
(u, T) — )■ (v, T) is equivalent to choosing the intersection of the asymptotes as the origin 
of coordinates. 

To show the existence of this scaling behavior, we compute the velocity without 
averaging over the disorder, i. e. we compute Vi for each of the = 14 samples that we 
have studied, with i = 1, ..,N. This choice allows to better see if all data scale on a single 
curve. Different system sizes = (100, 400, 1000)Al i.e. = (5000,8000,20000), 
different pinning strengths 10^ x ap/a^, ~ (2,3,5,7) and different realizations of the 
random position of the pinning centers have been investigated. In Fig. [TD] we plot 
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10"^ 10"^ 10"^ 10° 10^ 10^ 10^ lO'^ 



T\f\ 



-^6 



Figure 10. Scaling plot of versus T|/|-'^'' with /3 = 0.29 ± 0.03 and 

6^^ = 0.28 ± 0.05 for = 14 samples with different pinning strengths 10"^ x ap/a^ « 
(2, 3, 5, 7), system sizes = (100, 400, 1000)Ai and disorder realizations. Each sample 
value of /3 and is allowed to vary within the error bars. 



Vi\f\~'^ with respect to T|/|~^'' and a good collapse of all individual samples to a single 
curve with two branches is obtained. The values of used to compute / are sample 
specific. The two prefactors Vq and Tq are not universal and their value varies for each 
sample, but the scaling functions S± remain unchanged. We verified that averaging over 
disorder as previously defined, we obtain the same scaling functions. 

The collapse of our data implies that a true critical regime is observed with power 
law responses for the velocity with respect to the driving force and to the temperature. 
Moreover the critical depinning exponent /3 and the thermal exponent do not depend 
on the strength and the positions of the pins. Theses results suggest some degree of 
universality of our model for the elastic depinning of vortices. 



5. Discussion 



In our study we find a continuous elastic depinning transition for weak pinning strength. 
Both the depinning exponent /3 and the thermal exponent S, and the scaling function 
linking velocity v, temperature T and driving force / show some degree of universality 
with respect to disorder. A large variety of /3 values has been found in the literature as 
cited in section [H Our value is close to the value /3 ~ 1/3 measured for interfaces [60] in a 
space dimension d = 2. However, our result should rather be compared to other periodic 
systems. For example, our result is close to the value /3 = 0.35 recently measured 
in anisotropic periodic systems (stripe systems) ^3] for which the displacement field 
dimension isA^ = 2ina(i = 2 dimensional space. Much studies appear in the literature 
for the case = 1, in particular for CDWs where numerical simulations [39] in d = 2 
give /3 = 0.65 ± 0.05 and = 0.5 ± 0.1, whereas our results give very different values 
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/3 = 0.29 ± 0.03 and u = 1.04 ± 0.04. Also, we find u ^ 1 for the intrinsic correlation 
legnth exponent while z/ = 1/2 is expected for CDWs above threshold. Moreover, our 
results of section 13.21 suggest u = ups in contradiction with CDWs. Actually, the poor 
agreement between our results and CDWs may originate from the fact that = 1 for 
CDWs whereas A^ = 2 in our vortex simulation. 

We note that inserting z/ = 1 in the scaling relation v = 1/(2 — C) which originates 
from the statistical tilt symmetry [20] gives ( = 1. When inserting both values in 
the hyperscaling relation /3 = ^{z — () [IHl EDI [63] (where z and ( are respectively 
the dynamic and roughness critical exponents) together with z = 5/4 [65], we obtain 
/3 = 1/4 which is close to our result /3 = 0.29 ± 0.03 but may not be satisfactory since 
C = is expected to hold in the hyperscaling relation. 

Finally, going back to the connection between the depinning transition and 
absorbing-state phase transitions, we note that our result (3 = 0.29 ± 0.03 is close 
to the value found in d = 1 directed percolation or conserved directed percolation, and 
that z/ = 1.04 ± 0.04 is close to the value of z^_l in (i = 1 directed percolation. 

6. Conclusion 

In this paper we have numerically studied the 2D superconductor vortex dynamics in 
random media. A crossover from plastic dynamics dominated by disorder to elastic 
dynamics dominated by elasticity is found. We investigated the depinning transition in 
the elastic regime. Above the depinning threshold all the vortices depin together and 
flow in elastically coupled rough static channels. Three dynamical regimes are found: 
the first one is controlled by the finite-size of the simulation box leading to the so-called 
single-particle regime, the second one is the critical regime from which a power law 
response w ~ at T = can be extracted with /3 = 0.29 ± 0.03, while the third one 
is the recovery of the ohmic regime. A finite temperature analysis gives the thermal 
rounding exponent = 0.28 ± 0.05 defined by the power law response v ~ T^/^ at the 
threshold force. A finite-size scaling analysis at zero temperature suggests (in contrast 
with CDWs) the existence of a unique diverging length above the depinning threshold 
with an exponent u = 1.04±0.04, which controls the critical force distribution, the finite- 
size crossover force distribution and the intrinsic correlation length. A scaling law for the 
temperature and the force dependence of the velocity is found, consistent with a second 
order phase transition. Both critical exponents and the scaling function are independent 
of the disorder (strength and position of the pins) in the range of parameters we studied, 
indicating some degree of universality. Nevertheless the comparison with other similar 
periodic systems may suggest that large universality can not be found for 2D elastic 
depinning. We expect that such elastic depinning transition occurs in relatively weak 
pinning regimes of superconductor vortices. Within transport measurements where 
the voltage response to an applied current is equivalent to the velocity response to an 
applied force, the critical exponents /3 and 6 may be accessible from finite temperature 
measurements at different driving currents. 
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